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ABSTRACT 

Context. Oscillations in coronal loops are usually interpreted in terms of uncoupled magnetohydrodynamic (MHD) waves. Examples 
of these waves are standing transverse motions, interpreted as the kink MHD modes, and propagating slow modes, commonly reported 
at the loop footpoints. 

Aims. Here we study a simple system in which fast and slow MHD waves are coupled. The goal is to understand the fingerprints of 
the coupling when boundary conditions are imposed in the model. 

Methods. The reflection problem of a fast and slow MHD wave interacting with a rigid boundary, representing the line-tying effect of 
the photosphere, is analytically investigated. Both propagating and standing waves are analysed and the time-dependent problem of 
the excitation of these waves is considered. 

Results. An obliquely incident fast MHD wave on the photosphere inevitably generates a slow mode. The frequency of the generated 
slow mode at the photosphere is exactly the same as the frequency of the incident fast MHD mode, but its wavelength is much smaller, 
assuming that the sound speed is smaller than the Alfven speed. 

Conclusions. The main signatures of the generated slow wave are density fluctuations at the loop footpoints. We have derived a 
simple formula that relates the velocity amplitude of the transverse standing mode with the density enhancements at the footpoints 
due to the driven slow modes. Using these results it is shown that there are possible evidences in the observations of the coupling 
between these two modes. 

Key words. Magnetohydrodynamics (MHD) — Waves — Magnetic fields — Sun: atmosphere — Sun: oscillations 



X 



1. Introduction 

In coronal loops there is evidence of the presence of fast stand- 
ing magnetohydrodynamic (MHD) modes and slow (propa- 
gating and standing) MHD modes. Stand ing kink osci l lations 
were first reported us i ng TR ACE by fAschwan den et al.l d!999t) 
and iNakariakov et al. ( 1999|). Later, simila r observations were 
analys e d by, e.g.,[Sc hi'iivei' & Brownl (l2000l) : lAschwanden et al] 
' d2002l) : ISchriiver et al.l d2002l) . There is also clear evidence of 
the presence of pro pagating slow waves at loop footpoints (see 
iDe MoorteH l2009i for a review). These slow waves are most 
likely due to coupling between the underlying atmospheric lay- 
ers since the dominant periods tend to be around 5 min, sug- 
gesting a possible link with solar p-modes. Furthermore, obser- 
vations of standing slow modes in coronal loops, with periods 
largely above 5 min, have been also reported by a number of au- 
thors, using different instruments such as SOHO/SUMER (e.g., 
Wan g"etail l2003allbl [2007b . Yohkoh/BCS ( e.g.. iMariskal [2 005. 
2006) and more recently, Hinode/EIS (e.g., Erd elvi & TarovanL 



The theoretical interpretation of such a variety of oscilla- 
tions is done, in most of the cases, in terms of uncoupled MHD 
waves. On one hand, transverse oscillations are identified as fast 
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kink modes in the zero-/? approximation, meaning that in the lin- 
ear regime there is no longitudinal velocity along the magnetic 
tube. On the other hand, the reported slow modes are associ- 
ated to acoustic modes, ignoring the coupling with the transverse 
motions. Although these identifications are useful and simple, 
we have to bear in mind that under realistic conditions MHD 
modes may couple. In general, the effect of gas pressure, in- 
homogeneities or boundary conditions lead to the coupling be- 
tween fast and slow modes. As we shall see, a clear example of 
coupling comes from the line-tying condition, which amounts to 
considering the photosphere as providing a complete reflection 
of any coronal disturbance impinging from above. This is justi- 
fied in most of the cases due to the large difference in densities 
between the photosphere and corona but it obviously neglects 
the important role of the transition region. 

In this paper we study the simplest configuration in which 
fast and slow modes are mixed to understand the basics of 
mode coupling. We use the simple idea that a fast wave that is 
obliquely incident on a boundary will generate a slow wave. As 
far as we know, this issue has be e n addressed in s l ightly differ- 
ent co ntexts by, e.g. JSteinl ll971h : IVasquezl dl990h : loiiver et al.l 
dl992l) . A remarkable work about magnetohydrodynamic waves 
in coro nal flux tubes including th e line- tying effect was carried 
out by iGoedbloed & Halberstadtl d!994l) . In that work, it was 
clearly stated that pure fast or pure slow modes do not exist in 
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a line-tied coronal loop. Here we follow a different approach to 
analyse this problem. Our interest is in the effects of the applica- 
tion of line-tying conditions on fast and slow modes and in the 
possible observational fingerprints of this coupling. 

2. Model and Dispersion relation 

We consider first the simplest magnetic configuration where the 
magnetic field points in the z-direction, the plasma density and 
pressure are constant, and Va > c s (va is the Alfven speed and 
c s the sound speed). Since the medium is homogeneous, we con- 
sider perturbations that are proportional to e ! (&»+fe*+* z z) i meaning 
that waves propagate in the negative x and z-directions. In this 
configuration, the linearised MHD equations (see Appendix [A} 
lead to the well known dispersion relation for fast and slow 
MHD waves, 



{k 2 K + k 2 ) (c 2 s + vi) a? + k 2 (k 2 x + k 2 ) c 2 v| = 0. 



(1) 



In our problem it is more convenient to assume that ai and k x 
are known and solve for k z . The reason is that the frequency of 
an incoming wave remains constant in the reflection problem 
while the longitudinal wavenumber can change. The dispersion 
relation written as a biquadratic equation for k z is, 



2 2 



k\ - \{c\ + ^o?-c 2 v\k 2 x \k 2 z 
[a> 2 -k 2 x (c 2 s+ vl)]=0. 



(2) 



We denote by kp and ks the fast and the slow longitudinal 
wavenumbers that are solutions to the previous biquadratic equa- 
tion, i.e., 



k 2 - 
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2 2 
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(3) 
(4) 

(5) 
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These wavenumbers are associated to the same frequency co and 
same k x , and can be approximated, us ing the small plasma-yS as- 
sumption, by (e.g. Oliver et all 1 1992b 



i-2 ^ ^L- k 2 
K F ~ 1 K x 

V A 

k 2 ~ — 

ci 



l 2 
C 

V V A 



(8) 
(9) 



These approximations will turn out to be useful in the following 
sections. According to Eq. <j8j the fast wavenumber may become 
purely imaginary for a certain choice of the parameters, indicat- 
ing that the wave is evanescent. Hereafter, we will restrict our 
analysis to propagating waves. 

It can be seen from the linearised MHD equations that the 
velocity polarisation for fast and slow modes is given by the fol- 
lowing expression 



v 7 = 



uj- 



c s k x k z 
-k 2 c 



(10) 



where k z is either kp or k$ • Note that if we change the direction 
of propagation along the field, k z changes to -k z and the polari- 
sation also changes sign. In our configuration and in the regime 
that we are interested (va > c s ), since kp < ks, fast modes are 
characterised by long wavelengths and by v x > v z , while slow 
modes have short wavelengths and v x < v z . The slow modes 
have a larger compression since the wavelength is small in com- 
parison with the fast modes. 

3. Propagating waves: the reflection problem 

3. 1 . Fast MHD wave reflection 

We consider the most elementary problem to show the process 
of mode coupling due to boundary conditions. A fast MHD wave 
travels downwards (in the negative z-direction) and represents an 
incoming wave. This wave interacts with the photosphere (lo- 
cated at z = 0, where line-tying conditions are applied) and re- 
flects (now travelling upwards). The amplitude of the fast incom- 
ing wave is F\, while the amplitude of the reflected fast wave is 
Fr. To satisfy the boundary conditions a slow MHD wave, also 
moving upwards, must be generated at z = 0. The excited slow 
mode has an amplitude Sq. It is easy to write the velocity com- 
ponents using the polarisation of fast and slow waves, given by 
Eq. ([Tol l, and the proper sign of the longitudinal wavenumber, 

v — c J(,oit+k x x+kpz) i j7 J(ojt+k x x-kfiz) 
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(12) 



Now boundary conditions are imposed at z = 0, representing the 
location of the photosphere. We use line-tying conditions, i.e., 
V x (z = 0) = V z (z = 0) = 0. According to Eqs. CTQ and O the 
following conditions must be satisfied, 



Fi + F R -S, 



or - kic. 



2„2 



cjk x k s 

2„2 



F\ - F R + Sq- 



k c c 



= o, 



= 0. 



(13) 



The solution in terms of the amplitude of the incoming wave 
(which is an arbitrary parameter) is 



Fi 
Fi 



(oj 2 - k 2 s cl)k F + (oj 2 - klc 2 )k s 

2c 2 k x k s k F 

(co 2 — k^c 2 ^ kp - (co 2 — k 2 c 2 ^j ks 



(14) 
(15) 



The coefficients depend on u>, k x (same for the three waves) and 
on longitudinal wavenumbers, kp and ks, which can also be ex- 
pressed in terms of u>, k x by using Eqs. (O and ©. 

When there is no coupling between fast and slow modes (i.e., 
when k x = 0) we have that |Fr/Fi| = 1 because u> 2 = k^c 2 , 
while Sc/Fi = 0. Hence, there is a complete reflection of the 
incoming fast mode while the slow mode is absent. When there 
is coupling, the reflection coefficient for the slow mode is always 
different from zero, meaning that the incoming fast MHD wave 



Terradas et al.: Coupling between fast and slow MHD waves 



3 



/ . 00 




' c Jv A =0.Z 


■ \ 


0.99 




c,/v A =0.4 




U. JO 








0.97 








0.96 




cyv A =0.6 




0.95 








0.94 









0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



Fig.l 

of the 
speed. 



, Amplitude of reflected fast mode given by Eq. ( 1 14b as a function 
horizontal wavenumber for three different values of the sound 



inevitably generates a slow mode at the boundary. In Figure Q] 
the reflection coefficient of the fast wave, i.e., Fr/Fi, is plot- 
ted as a function of k x for different values of the sound speed. 
The coefficient is close to unity for values of sound speed low 
in comparison with the Alfven speed, meaning that reflection is 
almost total. However, for relatively large values of c s /va the 
curves clearly show a minimum indicating that the reflection of 
the fast wave is less efficient. 




Fig. 2. Amplitude of the generated slow mode given by Eq. d!5t as 
a function of the horizontal wavenumber for three different values of 
the sound speed. The dashed line represents the approximation given 
by Eq. dl8t while the dotted line shows the position of the maxima 
according to Eqs. l |19l > and ( 120b . 



In Figure |2] the reflection coefficient of the slow wave, i.e., 
Sg/F\, is plotted. The behaviour of the reflected slow mode 
amplitude is basically opposite of the reflected fast wave. The 
curves show a maximum at the locations where the generation 
of the slow mode is more efficient. Again the value at the max- 
imum strongly depends on the ratio of the sound speed to the 
Alfven speed, and corresponds to a k x which is very similar to 
kp. The location of the maximum can be analytically determined 
from Eq. (T5[ and using the approximation of the slow mode 



frequency, u 2 

S g 2k x kp 
~F~i 



k 2 c 2 



(16) 



Since the frequency is fixed, we have from the approximate ex- 



pressions ©-(O for k 2 and ki the following relation 



ci 



s J 



+ k 2 ^ 



Inserting this expression in Eq. ( [ToT l we find that 

Sc 2c 2 k x k F 



kl{cl-v 2 A )-k 2 x {cl + V \) 



(17) 



(18) 



This is a good approximation to the slow mode reflection coeffi- 
cient (see dashed line in Figure|2]i, at least for the case c s <sc va. 
It is used to determine the location of the maxima of the curves 
in Figure [2] Imposing that the derivative of Eq. dTSl is equal to 
zero it is found that the maximum corresponds to 




v 2 A + c 2 ' 



(19) 



and the value of the reflection coefficient for this horizontal 
wavenumber is 



Sg 
Fi 



(20) 



In Figure|2]the location and value of the maximum is represented 
by a dotted line. There is excellent agreement between the ap- 
proximation and the location of the maximum calculated using 
the original expression for the reflection coefficient of the slow 
mode (Eq. (fT5ll). From Eq. (l20b we see that the maximum am- 
plitude of the generated slow wave depends only on the /3 of the 
plasma. 

Once we know the coefficients it is straight forward to calcu- 
late magnitudes of interest such as the density perturbation of the 
slow reflected wave as a function of amplitude of the incoming 
fast wave. As we will discuss later, these magnitudes can be di- 
rectly related to real observations. From the continuity equation 
(see AppendixlAl the density perturbation for the slow mode is 



= - (k x V x + k z V z ) = — 
tjj a) 



1 - 



• k 2 r 2 ^ 



k§c s 



S G - 



Using again the fact that a? k$ c 2 and Eq. ( 120b we find 



(21) 



(22) 



This equation relates the maximum density perturbation associ- 
ated to the generated slow mode with the velocity amplitude of 
the incident fast wave. 



3.2. The time-dependent problem 

The time-dependent problem is numerically solved to demon- 
strate the mode coupling described in the previous section. The 
temporal evolution of the waves and specially their interaction 
with the boundaries provides a clear picture of the coupling and 
complements the analytical results. 
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In this numerical experiment a pulse in the v x component is 
generated at t — 0. This pulse has the following form 



t/r A = 3.03 



V X (Z, t ■ 

v z (z, t ■■ 



0) = vocos(fe) e ( « ) 
0) = 0, 



(23) 
(24) 



and the rest of the perturbed variables are set to zero. This 
particular profile has the property that represents a localised 
wave packet but with a rather well defined wavelength (basically 
given by A — 2n/k). To show the interaction of the fast mode 
with a single boundary, rigid conditions are applied at z — 
(v x (z = 0) = v z (z — 0) = 0) while open conditions are imposed 
at z = L. The linearised time-dependent MHD equations (see 
Appendix[A} are numerically solved using standard finite differ- 
ences techniques. 

The temporal evolution of the two velocity components is 
plotted in Figure [3] The initial disturbance excites the fast mode 
but also the slow mode since the initial perturbation does not 
satisfy the velocity polarisation relation, given by Eq. ( TTOb . The 
excited fast and slow modes split in two identical modes prop- 
agating in opposite directions (top panel). The fast and slow 
waves travelling to the left are denoted as Fj and Si, respec- 
tively, while F\+ and S\+ move to the right. Once the F\ mode 
reaches the rigid boundary (z = 0), it gets reflected (see mid- 
dle panel) with an amplitude Fr and a slow mode is generated 
(see dotted line), i.e., the Sc mode according to the notation in- 
troduced above. The slow wave packet has a typical wavelength 
much smaller than the fast wavelength {kp <s ks). At the other 
boundary (z = L) the fast mode F\+ leaves the system due to the 
transparent boundary conditions. At later times (bottom panel) 
the Sc and S\ modes, moving in opposite directions, start to su- 
perpose, while the and Si+, travelling in the same direction 
also interfere. 

The amplitudes of reflected fast and slow waves estimated 
from the time-dependent problem agree quite well with the cal- 
culations based on the analytical reflection coefficients given by 
Eqs. (fT4l-(fT5]l. Note that these equations were derived for purely 
sinusoidal and monochromatic waves but the initial packet is 
well represented by a dominant frequency (wavelength), and this 
is the reason of the agreement with the analytical results. 

3.3. Slow MHD wave reflection 

Although we are mainly interested in the fast reflection prob- 
lem, for completeness the slow mode reflection problem is also 
studied. The difference with respect to the fast MHD reflection 
problem is that reflection of the slow MHD mode at the bound- 
ary generates a fast MHD mode. The velocity components for 
this problem are 



/,|2 _ lA r 2 
V = Si e K<ot+k x x+k s z) 

c 2 k x k s 

+ Fq ^fi >tJ !-k x x—k^£) 



z) -s R 



/,,2 _ i-2 2 



cjk x k s 



y — g l e i(.<t>t+k x x-k$z) _|_ £ e i(u>t+k x x-ksz) 

c s k x kp 



Fq- 



i(b)t+k x x—kpz) 



(25) 



(26) 



Applying rigid boundary conditions the following coefficients 
for the reflected slow mode and the generated fast mode are 
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8 10 



t/T,= 7.58 




8 10 



t/r, = 12.12 




Fig. 3. Time evolution of the propagating disturbance given by 
Eqs. J23 M24t , The continuous line represents the v v component while 
the v z component is plotted with a dotted line (its amplitude has been 
multiplied by a factor 40 for visualisation purposes). Fast and slow 
modes travelling to the left are denoted as Fi and Si, Fj+ and S ;+ move 
to the right, and S G is the generated wave at z = (moving to the right). 
For this simulation zo/a = 5, ka = 27r/0.75, vo/v A = 1, k x a = 1.5 and 
c s /Va = 0.2. In this plot the time has been normalised to the character- 
istic time scale, t a = va/o.. 



(27) 



(28) 



found 

S R (a) 2 - k 2 s c 2 ) k F + [cj 2 - klc 2 ) k s 

Si (w 2 - klc 2 ) k F - (w 2 - klc 2 ) k s ' 

F c 2 (oJ 2 -k 2 cl){a) 2 -klcl)k s 

Si ~ c 2 A,{a) 2 -k 2 c 2 )k ¥ -{a) 2 -klc 2 )h 

It is worth to note that the reflection coefficient of the slow mode 
is exactly the same as the reflection coefficient of the fast mode 
given by Eq. (TT41) . 

4. Standing waves 

Once we understand the reflection of a propagating fast and a 
propagating slow wave at a rigid boundary we extend our anal- 
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ysis to a more realistic physical situation, i.e., the standing fast 
wave problem that is commonly reported by TRACE. 

There are several ways to study this problem, and the per- 
turbation method is one of them. Under this approach, we know 
that for zero-/? we have a pure standing fast wave, with only a 
v x velocity component, that satisfies the boundary conditions. 
When /3 is different from zero but small, a v z component is in- 
troduced, due to Eq. ( fTOb . and even worse, this component does 
not satisfy the boundary condition. In order to compensate, we 
have to add to the system a slow wave perturbation (with a dom- 
inant v z in comparison with v x ), which combined with the v z of 
the fast, will satisfy the boundary conditions jointly. Now the v x 
introduced by the slow mode is irrelevant since it is higher order 
in /3 (which is assumed to be small). 

Another way to analyse this problem is to solve the ful l 
stan ding problem. This was already d one bv lOliver et al 1 (119921) 
and iGoedbloed & H alberstadt ( 1994) and it is based on the su- 
perposition of the fast reflection problem plus the slow reflec- 
tion problem. The perturbation scheme and the full eigenmode 
problem give similar results, except where the solutions are 
completely mixed, i.e., when the slow mode component is not 
small in comparison with the fast component. This takes place 
around the "avoid ed crossings" in the dispersion diagram (see 
lOliver et al.L ll992l for further details). 

The advantage of the perturbation scheme is that we can 
still obtain simple approximations for the velocity and density 
changes around the footpoints. Let us assume that an initial 
disturbance mostly excites the standing pattern in the v x com- 
ponent. This standing wave is composed by the superposition 
of two identical propagating waves travelling in opposite direc- 
tions. We concentrate on the fundamental mode, having a maxi- 
mum at z — L/2 and a node of the velocity at the boundaries. A 
propagating slow wave is induced by the reflection at the bound- 
aries of the incident propagating fast wave that forms part of the 
standing fast pattern. Hence, if the transverse velocity amplitude 
of the standing fast wave at the loop apex is Vo the amplitude 
of the incoming fast wave is Vo/2 (the reflected fast wave will 
also have an amplitude of Vo/2 in our approximation). Once we 
know this amplitude it is easy to estimate the velocity or density 
fluctuations associated to the slow mode using the results of the 
propagation problem. According to the analysis performed in the 
previous section (see Eqs. d20l ) and (1221) ) 



V 



and 

Pi 
P 



(29) 



(30) 



These very simple expressions give the order of the maximum 
velocity and density perturbations associated to the slow re- 
flected mode as a function of the transverse velocity amplitude at 
the loop apex and the characteristic speeds of the configuration. 
When c s <sc va the velocity perturbation reduces to 



cl 



Vo 



1 



- 1 max j /-> a i 

v\2 4 



while the density perturbation is 



c s Vo 1 _Vo 

-T^r = -jyfi—- 

v. 2 4 c s 



(31) 



(32) 



It is worth to calculate the order of magnitude of velocity and 
density fluctuations around the footpoint. Let us take typical 
values: v A = 800 km s _1 , c s = 200 km s -1 , V = SOkms" 1 . 
Thus, the velocity perturbation according to Eq. ( 1291 is |v,| max =s 
2.5kms~'. Using Eq. ( l30l l we find that the density fluctuation 
associated to the slow mode is |pi/p| max — 1.25%. These den- 
sity perturbations, although small, should be detectable with the 
current instruments. Obviously, larger sound speeds and or small 
Alfven velocities would result in an increase of velocity and den- 
sity fluctuations. 

The previous estimations are based on the assumption that 
the excited slow modes do not have time to reach the opposite 
footpoint and reflect. If this is the case, the problem is more diffi- 
cult since the slow modes interfere and the velocity and density 
changes can be larger than those given by Eqs. (I29l>-(f30b. It is 
also important to remark that the density perturbation associated 
to the fast standing mode has the same profile as the velocity (the 
v x component), therefore the fast density perturbation is zero at 
the footpoints (and maximum at the apex), this explains why it 
is enough to consider only the density changes associated to the 
excited slow modes around the footpoints. 

4. 1 . The time-dependent problem 

Now the time-dependent problem is solved. The initial perturba- 
tion, representing the excitation of a mainly fast standing wave, 
has the following profile, 



v x (z, t = 0) = v sin (kz) , 
v z (z,t = 0) = 0. 

4.1.1. Rigid boundary conditions 



(33) 
(34) 



Line-tying conditions are applied at the two loop footpoints. We 
choose the longitudinal wavenumber that excites the fundamen- 
tal standing fast MHD mode, i.e., k = n/L. The results, repre- 
sented in Figure|4]manifest the generation of slow modes at the 
footpoints (top panel). The fast mode drives slow modes which 
move in the direction of the loop apex, located at z = L/2 (see 
middle panel). At a certain time (when t L/2c s ) the interfer- 
ence between the opposite propagating slow waves is produced. 
Since the two slow waves are identical but travelling in opposite 
directions a quasi-standing pattern is visible (see the profile of 
v z around z = 5a in the bottom panel). Eventually, slow modes 
reach the opposite footpoint and get reflected. Under such con- 
ditions the inverted process takes place, i.e., the slow mode gen- 
erates a fast mode. Nevertheless, this problem would take us too 
far and it is out of the scope of this work. However, it is worth 
to note that depending on the equilibrium and wave parameters 
the generated slow modes might match the frequency of a stand- 
ing slow eigenmode of the loop. This means that since the slow 
mode is driven by the fast standing mode, the frequencies of the 
slow standing and fast standing waves will be basically the same 
and the modes will have a highly mixed nature. This takes place 
around the "avoided crossings" in the dispersion diagram. 

4.1 .2. Sharp photosphere-corona transition 

From the previous results we might think that the generation of 
slow modes takes place only when rigid boundary conditions are 
imposed. In order to address this question, instead of line-tying 
conditions, a density transition between the corona and the pho- 
tosphere is considered. As density profile we use the following 
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Fig. 4. Time evolution of the standing fast excitation given by 
Eqs. J33M341. The continuous line represents the v v component while 
the v ; component is plotted with a dotted line (its amplitude has been 
multiplied by a factor 10 for visualisation purposes). For this simulation 
ka = n/10, vo/va = 1, k x a =1.5 and c s /va = 0.2. 

idealised model, 

P = (pph-p c )e~(») + (p p h -p c )e -( ^) +p c - (35) 

The parameter w measures the width of the transition between 
the photosphere, with density p p h, and corona, with density p c . 
The density contrast is set to p p h/p c = 10 8 . Non-reflecting condi- 
tions are imposed at the boundaries, now placed below the pho- 
tosphere. The details of the model used to represent the transi- 
tion photosphere-corona are not important since our focus is on 
the generated slow modes. The equilibrium gas pressure is set to 
the same constant value in the corona and photosphere to have 
magnetohydrostatic equilibrium. 

The results of the time-dependent simulation are plotted in 
Figure [5] There are some differences with respect to the be- 
haviour found for the perfectly reflecting boundaries. For exam- 
ple, the system now allows the energy to escape through the pho- 
tosphere, i.e., there is a transmitted fast (and slow wave). This 
means that the energy can leak through the photosphere but this 
is not visible in Figure [5] because the amplitude is very small. 
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Fig. 5. Time evolution of the standing fast excitation given by 
Eqs. l!33 M34t with the density model given by Eq. 1351 . The continuous 
line represents the v,- component while the v, component is plotted with 
a dotted line (its amplitude has been multiplied by a factor 10 for visu- 
alisation purposes). The vertical dashed lines represent the locations of 
the transition between the corona and photosphere. For this simulation 
ka = n/9, Vo/va = 1, k x a = 1.5, cJv A = 0.2, p ph /p c = 10 s , w/a = 0.25. 

Nevertheless, the primary result here is that a similar slow mode 
generation in comparison with the case of reflecting boundaries 
is produced in the coronal part (compare Figures [4] and 0. The 
results are qualitatively similar, thus, the overall conclusion is 
that nonuniformity causes the coupling between fast and slow 
modes in a similar way as line-tying conditions. 

5. Link with observations 

It is interesting to look for signatures of the mod e coupling stud- 
ied in this work in the observations. Recently, IVerwichte et al.l 
d2010h . have reported, using combined observations from 
TRACE and EIT/SoHO, intensity fluctuations at one loop foot- 
point with basically the same period, around 40 min, as the fast 
transverse oscillation of the loop. This long periodicity indicates 
that the origin of these density oscilla tions is most probably not 
photospheric. lVerwichte et all (l2010h suggest that the intensity 
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oscillations are due to variations in the line of sight column depth 
produced by the changes in the loop inclination as it oscillates 
with the transverse kink mode. This would explain the coinci- 
dence in periods, and also in the damping times. 

According to our study an alternative expl a nation for the 
intensity variations reported by IVerwichte et al.l d2010h is done 
in terms of the coupling between fast and slow waves due to 
the line-tying conditions. The similarity of the periods of the 
transverse mode and the intensity fluctuations is explained by 
our model where the excited slow wave, and the correspond- 
ing intensity oscilla tion, has the same pe riodicity as the trans- 
verse loop motion. IVerwichte et al.l (1201 Oh show intensity varia- 
tions associated with a transverse oscillation in a large loop of 
690±60 Mm length. The oscillation, which is a horizontally po- 
larised fundamental kink mode, has a period of P=2418±5 s. 
Therefore, the phase speed is equal to v p h=580±50 kms" 1 . For 
a thin loop in the long wavelength limit, this phase speed tends 
to the kink speed, Ck, of a cylindrical tube. If the intensity oscil- 
lation is associated with a slow magnetoacoustic mode with the 
same periodicity as the kink mode, then its wavelength is equal 
to A s =c s P=(c s P/v p h)2L. For a temperature between 0.9 and 1.8 
MK, the coronal sound speed is in the range 144-198 kms -1 . 
Hence, /l s =(0.30+0.07) 2L. For a standing mode, this translates 
into a wavenumber«=v p h/c s =2L//l s =3. 3+1.0. The EIT observa- 
tions suggests that the intensity variations at the two foot points 
of the loop oscillate in anti-phase. Hence, the mode, if stand- 
ing, should be an odd harmonic. Figure [6] shows the relative 
profile of intensity variations as a function of loop distance. 
The relative intensity has an amplitude of approximately 13%, 
which translates into a density perturbation with an amplitude, 
Ipi /plmax=6.5%. However, it is clear that a third or fifth harmonic 
could fit the observed intensity variation if we take into account 
that most likely the slow mode is still evolving towards a stand- 
ing wave pattern. We have measured the full intensity profile as 
a function of distance along the loop using EIT. Nevertheless, 
this has large uncertainties attached to it because of the lack of 
resolution and line-of-sight confusion. Also, the loop may be 
longitudinally structured in temperature. The EIT measuremen t 
hints at the presence of a fifth harmonic. IVerwichte et al.l {2010) 
show that the displacement amplitude, £o> can be modelled by a 
ten degree inclination of the loop. As the loop has a height of 236 
Mm, this means that £o=41 Mm. Hence, the velocity amplitude 
V()=£o<^=106 kms -1 . Using Eq. ( l30l l. and the estimated value for 
the Alfven speed, v A =410±40kms~ 1 , we find |pi/p| max =5+3%. 
This is consistent with the observed density amplitude (around 
6.5%). 

In another work. IVerwichte et al.l (12009) also reported inten- 
sity oscillations associated with a kink mode seen by EUV on 
board of STEREO and interpreted these as resulting from varia- 
tions in the line-of-sight column depth due to the loop showing 
a varying aspect to the observer. The measured relative density 
amplitude is 2.5%, and located around to loop top. Can this ob- 
servation be explained by a slow mode instead? In that obser- 
vation, the key parameters that were measured are L=340±15 
Mm, v ph =1100±100kms- 1 , V =35±9kms- 1 and v A =800+100 
kms" 1 . For the peak temperature of the EUV 171A bandpass of 
0.9 MK, c s = 144 km s~' . The slow mode would have a wavenum- 
ber n-l± 1 (in the case of a standing mode), corresponding to a 
wavelength of about /l s =90±10 Mm. Again, using Eq. (1301 . we 
calculate |pi/p| m ax=0.4%, which is negligibly small. Therefore, 
for this observation, the slow mode is not able to explain the 
observed intensity va riations. The main diffe rence with respect 
to the case studied bv IVerwichte et al] d2010h is that the density 



perturbations are located around the loop top instead of the loop 
footpoints. 
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Fig. 6. Relative profile of intensity variations se en in a transversely 
oscillating loop studied by Ver wichte et al] i2010h using TRACE as a 
function of distance along the loop. The solid circle is the measurement 
of the intensity variation at the loop foot point. The long-dashed and 
dashed curves are functions -0. 13sin(37r/L) and -0.13sin(57r/L), respec- 
tively, representing a possible third or fifth standing harmonic of the 
slow mode. 



6. Discussion and conclusions 

In this work we have investigated the possible effects of linear 
coupling between fast and slow modes that takes place due to the 
reflection at the photosphere. We have shown, by solving the re- 
flection problem, that the transverse motion of the loop produces 
slow MHD waves as long as the fast wave is obliquely incident 
on the boundary. These slow waves are manifested as propagat- 
ing density fluctuations with the same frequency as the standing 
transverse mode and are generated at the loop footpoints. The 
wavelength of these slow modes is basically given by k s ^ a)/c s , 
and since v\ > c s their wavelength is smaller than that of the 
corresponding fast standing transverse mode. Moreover, we have 
derived using the reflection coefficient of the slow mode, simple 
analytical expressions that relate the transverse displacement at 
the loop top with the amplitude of the density fluctuations pro- 
duced at the footpoints. The coupling between fast and slow is 
proportional to the plasma-p 1 , indicating that under coronal con- 
ditions it will be weak in general. However, it is proportional to 
the amplitude of oscillation of the fast mode, which can be large 
in some cases, and close to the footpoints, i.e., close to the pho- 
tosphere, the ratio between the sound and Alfven speeds can in- 
crease. Ideally, from the properties of the reflected slow modes, 
we could have indirect information about the real line-tying con- 
ditions in coronal loops and realistic values of the sound speed 
near the footpoints. 

The generated slow waves can eventually form a standing 
pattern along the loop. Thus, mode coupling might provide an 
excitation mechanism of standing slow modes, different from the 
driven photospheric origin usually emphasised in the literature. 
However, since the sound speed is assumed to be smaller than 
the Alfven speed, the time required for the slow mode to travel 
along the loop and reflect at the footpoint is much larger than that 
of the fast MHD mode. This means that slow standing modes 
will need much more time than fast standing modes to build up. 
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Nevertheless, if the wavelength of the slow modes is very short 
they might be damped by thermal conduction before the standing 
wave is formed. 

It has been show n that in the observations analysed by 
IVerwichte et al.l (120101) there are possible fingerprints of the cou- 
pling between fast and slow modes due to line-tying, providing 
an alternative explanation to the effects of integration along the 
line of sight pro posed by these authors. However, in the observa- 
tions studied by Verwic hte et al.l ((2009) the mechanism is unable 
to explain the intensity oscillations. One of the reasons might 
be that the density changes observed in this event are reported 
around the loop apex, while the estimations are based on prop- 
agating waves around the footpoints. In the case of the standing 
pattern the density fluctuations can be much larger. Further anal- 
ysis of other observations will help to test the linear coupling as 
an operative mechanism in coronal loops. 

It is worth to mention that the mode coupling discussed 
in the work is a purely linear effect. It is unrelated to the 
nonlinear coupling betwe en fast and slow modes due to the 
ponderomotive force (se e iHollwegl Il97lt iRankin et all 1 19941 
Terradas & Ofman, 2004). or th e parametric coupling studied by 
Zaqarashvili et al l d2002l 120051) . It is also different from mode 



conversion that takes place when c s = va. Line-tying boundary 
conditions are the responsible ingredients of the coupling, but 
we have shown that they are not the only way in which fast and 
slow waves may couple. A sharp transition between the corona 
and the photosphere produces the generation of slow modes, i.e., 
a change in the properties of the medium, for example the tran- 
sition region, leads inevitably to the coupling. 

Our theoretical model is based on a Cartesian geometry with- 
out a density enhancement representing a loop. For this reason, 
it is necessary to improve the model by including a slab or a 
cylinder that mimics the effect of a dense magnetic tube. This 
will complicate the mathematical problem of the reflection at 
the boundary, and it might be difficult to derive simple analyti- 
cal expressions. For example, in a cylindrical model, the eigen- 
functions of the slow and fast MHD waves have different radial 
dependencies and might share the same azimuthal wavenumber 
(playing basically the role of k x in the present work). This issue 
will be addressed in a future study. 
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Appendix A: 

We use the linearised MHD equations of the momentum, induc- 
tion, energy and continuity equation. Assuming Fourier analysis 
in the x— direction we have: 



dv x 
dt 

dvz 
~~dl 
db x 
~dt 

db, 

dt 
dpi 

dt 
dpi 

dt 



1 / ., ., B B db x 

- -iKxPi - ik x — b z + — 

P\ H P- oz 

_\_d_Pi_ 

p dz ' 
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-B ik x v x , 



JP ik x v x + 



dv z 
~dl 



~P I ik x v x + — 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 



In these equations Bo is the equilibrium magnetic field and p 
and p the equilibrium density and gas pressure, respectively. The 
rest of the variables correspond to the perturbed magnitudes. To 
eliminate the imaginary complex numbers we simply define 



v x = i v x , 
b* x = i b x . 



(A.l) 
(A.8) 



With this transformation the time-dependent equations are 
solved using standard numerical techniques. 

The dispersion relation, given by Eq. ([1]), is easily derived 
assuming in the previous equations a temporal dependence of 
the form e'°" and a spatial dependence with the z-coordinate of 
the form e' k ~ z . We have used the following standard definitions 
for the Alfven and sound speeds, 



V A = 



Cs = 



B 



(A.9) 
(A. 10) 



